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Abstract 

Analyses of observable supernova neutrino oscillation effects require the cal- 
culation of the electron (anti)neutrino survival probability Pee along a given 
supernova matter density profile. We propose a simple analytical prescription 
for Pee 5 based on a double-exponential form for the crossing probability and 
on the concept of maximum violation of adiabaticity. In the case of two-flavor 
transitions, the prescription is shown to reproduce accurately, in the whole 
neutrino oscillation parameter space, the results of exact numerical calcula- 
tions for generic (realistic or power-law) profiles. The analytical approach 
is then generalized to cover three-flavor transitions with (direct or inverse) 
mass spectrum hierarchy, and to incorporate Earth matter effects. Compact 
analytical expressions, explicitly showing the symmetry properties of Pee, are 
provided for practical calculations. 



PACS number (s): 14.60.Pq, 97.60.Bw 
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Typeset using REVT^ 



I. INTRODUCTION 



Observable effects of supernova neutrinos in underground detectors represent a subject 
of intense investigation in astroparticle pliysics, botli on general grounds (see tfie reviews in 
imj^) and in relation to the SN1987A event (see for an updated analysis and bibliography). 
In particular, flavor oscillations in supernovae may shed light on the problem of neutrino 
masses and mixing, by means of the (potentially strong) associated matter effects (see, 
e.g., |l|J5-0] for reviews of early works, and |]7|-p0[] for an incomplete list of recent studies). 
In particular, dramatic effects on oscillations have been predicted, related to the type of 
neutrino mass spectrum hierarchy and to Earth matter crossing. 

Given the importance of supernova neutrino oscillations for both particle physics and 
astrophysics, it would be desirable to have a simple and complete description of the most 
important quantity involved in the calculations, namely, of the electron (anti) neutrino sur- 
vival probability Pee-Q Approximate treatments of Pee have been proposed in the literature 
to cover parts the parameter space in a piecewise fashion, e.g., by using either the adiabatic 
approximation or the so-called resonance condition (see ||2l|J^, p!0| , p]T[] for recent examples), 
with inherent limitations in the range of applicability. In particular, it has been recently 
realized, first in the context of solar and then of supernova 0,0] neutrinos, that 



the time-honored resonance condition cannot be meaningfully extended at large neutrino 
mixing. Thus, apart from brute- force numerical calculations of Pee (see, e.g., p^p!8| -p0[]). 



a truly unified approach, valid in the whole three-flavor oscillation parameter space and 
applicable to generic supernova density profiles, seems still lacking in the literature, as far 
as we know. 

In this work, we propose a simple, unified analytical approach to the calculation of Pee, 
based on a double-exponential form for the crossing probability [^,^] and on the condition 
of maximum violation of adiabaticity [0,^ (replacing the popular resonance condition). In 
the case of two- flavor transitions (Sec. II), our prescription is shown to reproduce accurately 
the results of exact numerical calculations for generic (realistic or power-law) profiles, in 
the whole neutrino oscillation parameter space. The analytical approach is then generalized 
to cover three-flavor transitions with (direct or inverse) mass spectrum hierarchy (Sec. Ill), 
and to incorporate Earth matter effects (Sec. IV). Compact analytical expressions, explicitly 
showing the symmetry properties of Pee and useful for practical calculations, are summarized 
in the final section (Sec. V), to which we refer the impatient reader. 

II. TWO-FLAVOR TRANSITIONS 

In this section we discuss numerical and analytical calculations of the survival probability 
Pee for neutrinos and antineutrinos, assuming two-family mixing between Vf, and another 
active neutrino (z/^ = v^, or A simple analytical prescription will be shown to reproduce 
very accurately the exact numerical results for Pee. 



^The key role of Peeiy) and Pee^y)-, related to the practical undistinguishability of supernova muon 
and tau (anti)neutrinos, is neatly discussed in 
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A. 2u transitions: Notation 



In the case of two-family v^. Va oscillations (a = /i or r), we label the mass (m) 
eigenstates (z/i, z/2) so that vi is the lightest, 

Am^ = m\-m\> Q , (1) 

and parametrize the mixing matrix U as 



Uei Ue2 \ _ ( COS 9 sin 9 

Ual Ua2 / \ - sin 61 COS 61 



(2) 



where 9 G [0, 7r/2]. In vacuum, Ue —>■ Va oscillations can be described in terms of the 
pathlength (x) and of the oscillation wavenumber 

k = Am^/2E , (3) 

E being the neutrino energy. In matter, the z/g — *■ I'a dynamics also depends on the Ue — Va 
interaction potential difference P? 



V{x) = ^/2GfN,{x) , (4) 
where Nf,{x) is the electron density profile. In appropriate units, 

^^"^ -7.57xl0-i;(a:)4^, (5) 



eVVMeV ' ' ' g/cm^ 

where p{x) is the matter density and Yf,{x) is the electron/nucleon number fraction. 

In supernovae, p(x) [and thus V{x)] can be approximately described by a power law, 



p{x) oc X PS|. In the present work, power-law potentials are parametrized as 



X 



V{x) = V, j , (6) 

where n = 3 (unless otherwise stated), and distances are conventionally reported in units of 
the solar radius, Rq = 6.96 x 10^ m. 

Figure 1 shows an example of "realistic" neutrino potential profile V{x) (dashed curve), 
as graphically reduced from the supernova simulation published in for the case of M = 
14.6 Mq, M being the mass of the ejecta. In the same figure, the solid line represents the 
best-fit power-law potential (n = 3), corresponding to take Vq = 1.5 x 10~^ eV^/MeV in 
Eq. (|^). For definiteness, we will use the realistic or the power-law curves in Fig. 1 as 
representative V{x) profiles for our calculations. However, our main results are applicable 
to generic supernova density profiles. 
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B. 2zv transitions: Neutrinos 



Following the calculation of Pee{i^) from the initial v state in matter to the final v 
detection in vacuum^ can be factorized as 

/ N /cos^e sin^^A /1-P,(z/) P,{u) \ /cos2^„ sin2^^„\ /l\ 

where Pc{^) is the so-called crossing probabihty for neutrino^ [Pci,^-') = P{t^2rn J^i)], and 
6m is the effective mixing angle in matter at the origin, defined by 

sin 2^ 

sin2g^ = ^ , 8 

(cos20-y/A;)2 + sin2 29 



cos29-V/k 

cos 29m = I = ■ (9) 

^(cos20-\//A;)2 + sin2 29 

Note that in Eq. it is understood that oscillating terms are averaged out, thus providing 
an incoherent u state at detection. 

The high supernova core density (at the start of neutrino free streaming) implies V/k ^ 1 
in Eqs. (H) and (y), so that sin26'm — and cos26'm — —1- From Eq. (|^, one can then 
reduce the calculation of Pge to that of Pc, 

Pee{y) = cos^ 9 P,{u) + sin^ 9 [1 - P,(z/)] (10) 
= U',,P,{u) + U!,[l~P,{u)] . (11) 



1. Numerical approach 

One possible approach to the calculation of Pc (and Pee) is the numerical integration of 
the neutrino evolution equations along the supernova potential profile, as advocated in some 
recent works [^^-|5D[, as well as in a few earlier ones (see, e.g., [2^|). For the purposes of our 



work, we have performed a numerical (Runge-Kutta) calculation for Pc, assuming the two 
potential profiles in Fig. 1. 

Figures 2 and 3 show our numerical results as dotted isolines for Pc(z^) in the mass-mixing 
plane {Am'^/E, tan^^), for the case of power-law and realistic potential, respectively.^ The 
"bumpy" structure of the realistic V{x) profile is reflected by the "wiggling" behavior of 
Pc in Fig. 3, leading to significant differences with the Pc isolines of Fig. 2 in part of the 



^The discussion of possible Earth matter effects is postponed to Sec. IV. 

^The cases Pc — and Pc 7^ discriminate adiabatic and nonadiabatic transitions in matter. 

^The solid curves in Figs. 2 and 3 correspond to our analytical approximation for Pc, as discussed 
later. 
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parameter space. However, one can note that, for ls.rn? /E 0, the detailed structure of 
V{x) is irrelevant, and Pc{i^) cos^ 9 in both Figs. 2 and 3 (extremely nonadiabatic limit). 

Although brute-force numerical calculations of Pc can provide, in principle, "exact" re- 
sults, it should be stressed that they do not represent an optimal and efficient approach 
in the case of supernovae. Computer routines for integration are tipically time-consuming, 
being required to track a large number of oscillation cycles along a potential profile spanning 
many orders of magnitude. Instabilities and inaccuracies in the numerical results can easily 
emerge for realistic potentials (as experienced by ourselves), as a consequence of sudden 
variations in the lowest order V{x) derivatives.^ Moreover, numerical integration produces 
additional but useless (unobservable) information on oscillating factors and phases, and 
is thus inefficient for practical purposes.^ Last, but not least, the uncertainties affecting 
simulated supernova density profiles make it preferable to perform several approximated 
(but quick) calculations of Pc for different trial functions V"(x), rather than a single (but 
time-consuming) exact numerical calculation for a fixed V{x). 



2. Analytical approach 



The previous discussion indicates that a handy analytical approximation to the numerical 
results for Pc, applicable in the whole parameter space, is highly desirable. We propose (and 
motivate below) the following analytical recipe in three steps: 
(i) Identify the point Xp where the potential equals the wavenumber, 



V{xp) = k ; 

(ii) calculate the so-called density scale factor r at that point. 



(12) 



r = — 



1 dV{x) 
V{x) dx 



n -1 



(13) 



{Hi) Insert the above r in the double- exponential parametrization for Pc{v) (originally derived 
in the context of solar neutrinos |25|J2F 



Priv) 



exp{27Trk cos^ 



1 



exp(27rrA;) 



(14) 



The results of such an exceedingly simple analytical recipe are shown as solid isolines 
for Pc in Fig. 2 (power-law case) and Fig. 3 (realistic case), to be compared with the corre- 
sponding dashed isolines (exact numerical results). The agreement between numerical and 



^Indeed, we think that some numerical artifacts (fake wiggles) might be present in the numerical 
calculations of Pee as graphically reported in 



^For instance, the authors of |18,19| need to time-average their numerical probabilities, in order 
to force an incoherent initial state for the implementation of Earth matter effects. 
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analytical estimates of Pc is extremely good in the whole parameter space, the difference 
being SPc = 2 x 10~^ in the worst cases.Q The final accuracy for Pge is even better, since 
Eq. (|10|) implies SP^e = \ cos2e\6Pc < 5P^. 

Notice that, in the exact power-law case [Eq. (||)], the calculation of r through Eqs. (|12D 
and ( |I3|) is trivial. 

In the case of realistic potential profile, the only modest complication is the numerical 
solution of Eq. ([T^) and the evaluation of the derivative in Eq. ([T3| ) for the given (tabulated 
or parametrized) function V{x). 

Our effective analytical prescription for Pc{i^) in supernovae [Eqs. (p!2D -(|l4D] stems from 
several recent improvements in the understanding of nonadiabatic transitions, as discussed 
below. Although such improvements have been mainly tested in solar neutrino oscillations, 
they are often applicable also to supernova u oscillations. 



3. Discussion of the analytical prescription 

Equation (0), originally derived for the solar neutrino (exponential) density profile 
|26[| within the unnecessary restriction 6 < 7r/4, was explicitly shown in pO| , |3T| to hold 
for 6 > 71 /A as well, especially for appropriately chosen density scale factors r [EH]. The 
double-exponential form of Pc for 6 > tt/A has also been recently applied to the transitions 
of high-energy neutrinos from the decay of hypothetical massive particles trapped in the Sun 
32 1, and to the transitions of supernova neutrinos [15-0. Such parametrization for Pc has 



thus several desirable properties: (i) It is a good ansatz (as originally advocated in |Q) for 
generic density profiles; (ii) It holds in both octants of 9; {in) It reproduces the extremely 
nonadiabatic limit at small k; and (iv) It reproduces the single-exponential, Landau-Zener 
(LZ) limit at small 9 (see, e.g., ||^). As a final remark on Eq. (0), it should be noted that 
the single- and double-exponential forms for Pc involve, in general, a function F{9) depend- 
ing on the potential profile (see Table I in [§^). Our choice in Eq. (0) corresponds to take 
F{9) = 1 — tan^ 9, which is the exact result for a solar-like (exponential) profile, and repre- 
sents the leading prefactor of F in a 1/n expansion^ for a supernova- like (power- law) profile 
33| . We have checked that the inclusion of the full (more complicated) expression for F{9) 



in the power-law case [^,§,|T^ does not lead to a significant improvement in the (already 
high) accuracy of Pc of our analytical prescription, and thus we advocate the simpler form 
for Pc given in Eq. (|l^) also for supernovae. 

The concept of a "running" density scale factor r = r{xp) [as in Eq. (|T^)] was also 
originally introduced in the context of solar neutrinos ||26[ , typically by calculating r at the 



^We have investigated a variety of supernova density profiles available in the published or unpub- 
lished literature, and obtained similarly good results (not shown). 



The exponential profile case can be seen as the n —>■ oo limit of the power-law profile case |33|. 
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point Xp 
e.g. 



Xres defined by tlie so-called "resonance" condition V^Xr 



Am^ cos 26 fsee, 



33 



Such a choice for Xp, although successful at relatively small 9, is clearly not 
applicable for 6 > it/A [0], and fails to describe correctly nonadiabatic transitions at small 
k, where Pc ^ aX 6 ^ tt/4 [^,^. For large 6, the resonance condition can be misleading, if 
not meaningless, and it is more appropriate to characterize Pc through the point xmva where 
maximum violation of adiabaticity (MVA) is attained ||2^ , [23| , p!6[| . Indeed, in the context of 
solar neutrinos, the prescription Xp = xmva for r{xp) is more accurate and physically more 
consistent than Xp = x^-es l2^- In the context of supernova neutrino oscillations, it has also 
been recently suggested that xmva niight play an important role as well P^JT?]] , although 
the authors of p6|JT7[] do not use the prescription r = r(xMVA); but make an improved WKB 
calculation of the LZ exponent for Pc (involving a numerical integration in the complex 
plane). However, we have verified that the prescription r = r(xMVA) gives very accurate 
results for P^ in the whole mass-mixing plane for supernovae (very close to the solid isolines 
in Figs. 2 and 3), making it unnecessary, in practice, to resort to WKB- inspired or other 
relatively complicated approaches. In fact, our Eq. (|1^) is just a suitable approximation of 
the MVA condition, as we now discuss. 

For a monotonic V{x) profile, the MVA point is uniquely defined in terms of the flex 
point of cos29m{x) 



' d? cos 26m{,x) 



dx'^ 



. 



/ x=xmva 

For a power-law profile as in Eq. (Bl), the above equation implies that 



^(a^MVA) = k ■ g{n, 9) 



where 



, cos26' 2 — n 

' ^ 2 l-2n \ 



1 



n 



cos 26 2 — n 



1 - 2n 



l-2n. 



(16) 



(17) 



(18) 



and the ib sign must be chosen so that > 0. For n close to 3 (supernova case), by keeping 
only the leading term in the expansion of the square root in g{n, 6), one obtains: 



2 
3 
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g^l + 
n = 4 — > 5' — 1 



cos 29 - 1 
10 

cos 29 - 1 



7 



(19) 
(20) 

(21) 



namely, (7 ~ 1 for n ^ 3 ± 1, up to ~ 10% errors. 

A fractional error 6 in the evaluation of V^(xmva) leads, for power-law potentials, to a 
fractional error S/n in the evaluation of r(xMVA)- Therefore, by setting in any case g = 1 
[and thus = A;, as in Eq. ([I2|) , rather than V = kg, as in Eq. {^i7\) ] we expect a mere 
few % variation in the running value of r, which is of little relevance when r is inserted in 
Pc through Eq. (0), as we have numerically checked. 
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In conclusion, the condition V{xp) = k in Eq. (|^) represents a good approximation to 
the MVA condition for supernova neutrino oscillations, and can be used to replace the time- 
honored (but inapplicable at large 6) resonance prescription.^ One has then to calculate 
[through Eq. ([I3|)] the corresponding (running) value of the density scale factor r = r{xp) 
[provided by Eq. ( [T5| ) for an exact power-law profile], insert r in the double-exponential 
form for Pc [Eq. and finally get Pee through Eq. (|10D- Our analytical prescription is 

applicable to generic (power-law or realistic) V{x) profiles in the whole 2z/ parameter space, 
with a typical percent accuracy in Pge- 

Figure 4 explicitly reports the results of such a prescription for Pge for the two profiles 
in Fig. 1, and for a representative supernova u energy [E = 15 MeV). Notice that the Pee 
isolines for the power-law profile (dotted curves) are simply calculated through elementary 
functions [Eqs. (|T2|)-(|1^)]. Notice also that the Pee isolines for the realistic and power-law 
cases in Fig. 4 appear to differ significantly in two regions: (z) at relatively small mixing 
(tan^ ^ ^ 0.1); and (ii) at nearly maximal mixing (tan^ ^ ^ 1) with Am^ ~ 10~^^^ eV^. 
In connection with solar neutrinos, this fact implies that the detailed supernova density 
profile can be relevant for the oscillation parameters corresponding to the so-called small 
mixing angle (SMA), low Am^ (LOW), and quasivacuum oscillation (QVO) solutions to 
the solar neutrino problem [^. Conversely, for oscillation parameters within the so-called 
large mixing angle (LMA) solution at Am^ ^ 10~^ eV^, or within the multiple vacuum 
oscillation (VO) solutions at Am^ ^ 10^^ eV^ |^, transitions in supernova matter become. 



respectively, purely adiabatic or extremely nonadiabatic, with no significant dependence on 
the details of the V{x) profile. 

A final technical remark is in order. For solar neutrinos, it was shown in [^ that the 



MVA-inspired recipe for P^ [r = r(xMVA)] has to be matched and replaced, at small k, with 
r = const [ pT| , P^ . The constant (limiting) value for r can be elegantly derived, through a 
perturbative approach in terms of an integral involving V{x) in the convective zone of 
the Sun, where V{x) experiences a sudden drop. Conversely, in supernovae, V{x) vanishes 
in a smoother way at large x, and the small-A; correction (r = const) becomes unnecessary 
in practice. We have numerically checked that such correction does not appreciably improve 
the accuracy of the simple MVA-inspired recipe in Eqs. ([T^)-(p^. Therefore, concerning the 
calculation of Pc in supernova neutrino transitions, we advocate the use of Eqs. (|12D-(|I4[) in 
the whole 2i/ parameter space. 



C. 2z/ transitions: Antineutrinos 

The extension of our analytical prescription from the neutrino case to the antineutrino 
case can be obtained through the replacement V/k ^ —V/k in Pee. By conventionally 
keeping V > (and 9 unaltered), this imphes 



^It is curious to note that, in the context of solar neutrinos, the relevance of the point where V = k 
(as opposed to the resonance point) was suggested in [35] and then abandoned in the literature. 
See also [22] for an updated discussion. 
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Peeiyl + Am^) = Pee{y\ - Am^) . (22) 

The change of sign of Am^ for neutrinos is equivalent to a swap of the mass labels 1 2, 

- Am^) = P,,{v\ + Am%^2 , (23) 

corresponding to Ul^ ^ U^2 sin^ 9 ^ cos^ 9.^ 

The above two equations imply that, for fixed Am^ > [Eq. (|^)], Pee(^) can be obtained 
from Eqs. (p!oD and (0) through a 1 ^ 2 swap, namely, 

= sin^ 9 Pc{v) + cos^ 9 [I- P^{v)] , (24) 

where 

P,{V) = Pc{v)l^2 

exp(27rrA; sin^ 6*) — 1 , . 

exp(27rrA;) - 1 ' ^ ' 

with r [Eq. ([T3|) ] to be evaluated at the same (6'-independent) point Xp defined for the 
neutrino case in Eq. (p!2|). Isolines of Pc{p) from Eq. ( pSj) are just the mirror images (around 
the axis tan^ 6' = 1) of those obtained for Pc(^) in Figs. 2 and 3. 

D. transitions: Summary 

A neat summary of the previous analytical results for z/ and V can be obtained by 
introducing a new notation, 

p± ^ exp(±27rrfccos^g) - 1 
" ~ exp(±27rrA;) - 1 ^ ' 

_/ Pc{^) (+)> 

~V-Pc{v) (-). ^'^^ 
In terms of P^, the expressions of Pee{i-') and Pee(i^) for 2v transitions are unified as 

Pie = U!, Pt + U^2 (1 - P^) , (28) 

where the + sign applies to neutrinos, while the — sign to antineutrinos. 

In the above equation, for any sign (±), both V and k are kept > 0. All physical 
cases are covered by letting 9 span its whole range [0, vr/2]. The density scale factor r 
[Eq. (|1^)] is calculated at the point Xp where V{xp) = k [Eq. (p!2[)]. The accuracy of such 
analytical prescription, as compared with exact numerical results, is at the percent level 
for both realistic and power-law density profiles, as demonstrated in Figs. 2 and 3 for the 
representative cases shown in Fig. 1. In the power-law case [Eq. (^)], the calculation of r is 
further simplified [Eq. ([I5|) ]. 



^'^Such well-known symmetry properties are repeated here, in preparation of the more compH- 
cated case. 
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III. THREE-FLAVOR TRANSITIONS WITH HIERARCHICAL MASS SPECTRA 



Supernova neutrinos can provide peculiar tests of three-flavor oscillations in matter, ow- 
ing to the wide dynamical range of V{x) in collapsing stars (see |]7|- [14| , p!8| -|20[| and references 
therein for recent 3z/ studies). It is thus important to generalize the previous 2i/ results to 
the case of 3z/ transitions, as we do in this section for the phenomenologically interesting 
cases with hierarchical mass spectra. We think it useful to review the derivation of P^^ (re- 
covering some known results) with no reference to the often (mis)used concept of "resonant 
transition." 



A. transitions: Notation and phenomenological input 

We assume mixing among three active neutrinos {u^, v^, i^r) and three mass eigenstates 
(z/i, z/2, z/3) through a unitary matrix]^ U . The matrix elements Uei relevant for P^^ are 
parametrized in terms of two mixing angles (0, u) = (6*13, 6*12) [|], 

U^^ = COS^ COS^ LU , (29) 

U^2 = cos^ sin^ uj , (30) 
U',, = sin^ . (31) 

The kinematical parameters are completed by two independent squared mass differences 
(5m^, m^). The dynamics is fixed by V{x), and the full supernova 3z/ parameter space S^i, 
is 

= (5m2, m^, LU, 0, V) . (32) 



Solar and reactor neutrino oscillation analyses suggest ||36|j37 



Sm'^ = \ml-ml\<7x 10"^ eV^ , (33) 
while atmospheric neutrino analyses indicate 

~ \ml - 2I ~ 3 x 10"^ eV^ , (34) 
thus favoring the so-called hierarchical hypothesis 

Sm^ -C , (35) 

very often used in the literature. 

Under the assumption of Eq. (|35|), we parametrize the mass spectra (up to an overall 
mass scale) as 



In the context of supernova neutrinos, U can be taken real without loss of generality. 
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+m2) , direct hierarchy , 



(mi,m2,m3) = <^ (36) 



5m I (5m 



2 ' ' 2 ' 



+ ^~) — m ) , inverse hierarchy 



where, conventionally, mg — > (< 0) identifies the so-called case of direct (inverse) 
hierarchy, while m^— > for both hierarchies. As far as 0, G [0, '/r/2], such a convention 
can be shown to cover all physical cases, both in vacuum and in matter pO|| . 

Besides Eqs. (^3|)-(^), a further phenomenological input comes from the combination 



of reactor and atmospheric neutrino data, providing 1 37,38] 



sin^ (j) = [7^3 < few % . (37) 

As a consequence of the hierarchical assumption of Eq. ( pSj) [and, to some extent, also 
of Eq. (|37|)], the 3z/ dynamics approximately reduces to the dynamics in two 2z/ subsystems, 
dominated by relatively low (L) and high {H) values of the matter density, according to 
parameter space factorization 

Ssu ^ ® = (5m^ uj, V cos^ (p) (m^ 0, V) (38) 

(see 1^ and references therein). The corresponding "low" and "high" neutrino oscillation 
wavenumbers are defined by 

kL = Sm'^/2E (39) 

and 

kn = m'^/2E , (40) 

respectively. 

B. Sv transitions: Neutrinos with direct hierarchy 

For neutrinos with direct mass hierarchy, the factorization of dynamics in Eq. (RSf) leads 



to 





PM = ( U!, , , Ul ) I Pl{v) ' 1 - Pl{v) 




/I 

X 1 - Ph{v) Ph{v) I I f/e\m I , (41) 
VO Ph{v) 

where Pl and Ph are the crossing probabilities for the low and high density transitions, 
respectively, and the elements t/gj ^j in matter are defined in analogy to Eqs. (P^)-(PT|), but 
with neutrino mixing angles a;^ and 0^, in matter given by 
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sin 2ujr, 



cos 2uJn 



and 



sin 20n 



cos . 



sin luj 

'cos luj — cos^ V I + sin^ 2u; 

cos 20; — l^/Z^L 
^cos 2uj — cos^ l^/A;^)^ + sin^ 2uj 



sin2(/) 

cos2u; - V/kny + sin^ : 

cos 20 — V/ kn 
'cos2uj -V/kny + sm'^ : 



(42) 
(43) 



(44) 
(45) 



at zeroth order in 5m? /m?' The (high density) initial condition V/ki^u ^ 1 leads then 
to cos 20m — —1 — cos2c(Jm and to the known expression 



Pee(^) = f/fl Pl{v) Ph{v) + Ul, [1 - Pl{v)] Ph{v) + Uj, [1 - Ph(^)] • 



(46) 



On the basis of Eq. (^) and of our 2z/ prescription in Eqs. ([T^ ) - ([I^) , PlW) and Ph{j^) 
can be analytically expressed as 



PlM 



and 



Ph{i^) 



exp(27rrifc2, cos^ cj) — 1 
exp(27rri/cL) - 1 

exp{2iTrHkH cos^ 0) — 1 
exp{27irHkH) — 1 



(47) 



(48) 



where 



1 dV{x) 



(49) 



V(x) rfx 

and the points xl and xh are defined by the V = k condition]^ [Eq. (0)] 

^(a;L,/f) = kL,H ■ (50) 



Equations (|46D-(pOD allow the calculation of Peeii^) in the whole 3z/ parameter space for 
direct hierarchy. 



^^In principle, Eq. ( p8| ) implies an effective potential V{x) cos^ for the L subsystem, and thus 
V{xl) cos^ = k^. However, we have checked that, within the phenomenological bound in Eq. (|37|), 
the resulting difference in the calculation of Pee is completely negligible. We prefer then the simpler 
condition V{xl) = ki, analogous to V{xh) = kn for the H subsystem. 



12 



C. 3i/ transitions: Antineutrinos with direct hierarchy 

The antineutrino case can be obtained in analogy to the neutrino case [Eqs. (^TD-( 
through the replacement V/ki^n —V/k^ H- In particular, the antineutrino mixing angles 
in matter are given by 

_ sin2c(j 
sin2tt;m = / , (51) 

(cos 2uj + V cos^ (p/k^Y + sin^ 2uj 



o- cos 2uj + V/kL , . 

cos ZuJm = I 1 (5zj 

(cos 2uj + V cos^ 0/ k^)"^ + sin^ 2u} 



and analogously for 0^. The initial high-density condition gives cos 20^ ~ +1 ~ cos2a;TO, 
and leads to 

Pee(l7) = Ul, [1 - Pl{V)] + Pl{V) . (53) 

Our analytical prescription for Pl{P) is then obtained, mutatis mutandis, from Eq. (p5|), 

„ ._^ exp(27rrLA;Lsin2^) - 1 

^'<''' = exp(2.r,t,) - 1 • C^^' 
with ri defined as in Eqs. (^9]) and ([50|). 



D. transitions: Neutrinos with inverse hierarchy 

For fixed mass spectrum and mixing angles, neutrino and antineutrino probabilities can 
be transformed one into the other by flipping the signs of V/ki^n or, equivalently, by flipping 
the signs attached to (5m^ and m^, while keeping V > 0. The ±5m^ sign flip is equivalent to 
the 1 <-> 2 swap of mass labels (as discussed for the 2z/ case), leading to U^i U^2 [namely, 
sin^ti; ^ cos^a;, with no change in or U^^]. The ±m^ sign flip is instead equivalent to 
swap hierarchy [see Eq. (pB])]. 

From such symmetry properties one obtains 

Peeil^ I + Sm^, —m^) = Pee{l^\ — , +m^) = Peeij' \ + , +m^)i„2 , (55) 

namely, the Ve survival probability for inverse hierarchy equals the Ve survival probability 
for direct hierarchy [Eq. (|5BD], under the substitution Ul^ ^ f/Jji 

Pee(^) = Uj, [1 - PL{y)] + Ul PL{y) , (56) 

with Pl{^) defined as in Eq. (^71) . 
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E. 3iy transitions: Antineutrinos with inverse hierarchy 



In analogy with the previous subsection, it can be easily realized that 

Peeij'l + Sm^, —m^) = Peeil^l — Sni^ , +m^) = Pee{l^ \ + Sui^ , +m^)i^2 , (57) 

which, applied to Eq. (^6l), gives 

PUV) = Ul^ Pl{v) Pniu) + f/fi [1 - Piin PH{y) + [1 - PH{y)] , (58) 



with Pl{i^) and Ph{i^) defined as in Eqs (|5J) and (^8|) , respectively. 

Notice that we have written Ph{^) and not Pnip) in Eq. (p8D, since the 1^2 swap 
makes no change in Ph, defined in terms of mg — ml 2 and of U^^. One can then drop the 
argument of Ph and simply writeQ 

Ph{v) = Ph{i^) = Ph . (59) 



F. 2>v transitions: Summary I 

The 3z/ results in Eqs. (ffBI), (pBf), (pB]), and (BH) can be summarized as 



P!: = UlX, + Ul^X, + Ul,X, , (60) 

where the coefficients Xi are given in Table I, in terms of Pl{^), Pl^J'), and Ph- 

The parametrization of P^^ given in Eq. (^0]) and Table I agrees with the results pre- 
viously obtained in [^0. In addition, however, we have provided explicit analytical ap- 
proximations for Pl [Eqs. (^7]) and (|5^)] and for Ph [Eq. (^S])], allowing a straightforward 
calculation of Pg^g'' in the whole mass-mixing parameter space, for generic (realistic or power- 
law) potential profiles. 

G. 2>u transitions: Summary II 

An alternative summary for Pg^g^ can be obtained by introducing, in analogy with the 2u 
case, the notation 

p± _ exp{±2nrLkL cos^ u) - 1 
^ ~ exp{±27rrLkL) - 1 ^ ^ 

-I (62) 



^^A formal distinction between Ph{i^) and Ph{t^) was kept in the notation of [^ and then dropped 



m 
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and 



p± ^ exp(±27rrj^fcHCOS^0) - 1 
^ ~ exp(±27rrHA;iT) - 1 ^ ' 

^ I exp(-27rrHfcH sin^ 0) (+) , ^^^^ 



where, in the last equation, we have used the phenomenological inputs in Eqs. (|34D and 
(^), implying that 2'nrHkH cos^ 1 for typical supernova potential profiles and neutrino 
energies. 

Equation (^) can then be written in the equivalent form 

= Uli Pl Ph + Ul^ (1 - Pt) Ph + (1 - ^1) , (65) 

where the sign assignment for P^u is given in Table II for the four possible combinations of 
neutrino types (z/ or v) and mass hierarchy (direct or inverse). 

The neat summary given in Eqs. (^), (|63|) , ( |65l) and in Table II makes the symmetry 
properties of P^^ rather transparent. Passing from u to 17, or from direct to inverse hierarchy, 
appears to simply require appropriate sign flips. 



H. 3i/ transitions: Representative Pee calculations 

Figure 5 shows a representative analytical calculation of P^^ in the slice {5m?, tan^u;) 
of the ?)V parameter space, which is relevant for the L transition in supernovae [Eq. (PB[)] 
and for its connection with solar neutrino oscillations. Calculations (solid lines)0 are made 
at fixed E = 15 MeV, = 3 x lO^^ eV^, and tan^ = 4 x 10"^ (P^ = 0.48), for the 
power-law profile in Fig. 1. The upper panels refer to v (left) and V (right) in the case 
of direct hierarchy. The lower panels refer to v (left) and V (right) in the case of inverse 
hierarchy. Notice that all calculations for this figure involve only elementary functions, the 
density scale factor r being explicitly given by Eq. ([T5| ) for a power-law profile. 

A glance at Fig. 5 shows two apparent symmetries: Pee{i-') and Peeij') for direct hierarchy 
look like the mirror image of Pee{p) and Pee('^) for inverse hierarchy, respectively. This 
(approximate) symmetry originates from the small value of U^^ used in the calculations of 
Fig. 5. Indeed, neglecting the third term (proportional to U^^) in Eq. (|65|), and using Eq. (pSf ) 
(with the identification P^ = P^), it is 

P!: ^ P| P^, , (66) 

namely, the 3z/ probability is obtained from the 2z/ probability through a modulation factor 
P^. Using Eq. (^) and the symmetry properties in Eqs. (|55| ) and (p7D, it follows that 

P!:{y, direct) ^ P+ P^,:{v) = P+ P!:{v),^2 ^ P!:iv, inverse)i^2 (67) 



The dotted curves in Fig. 5 include Earth matter effects, to be discussed in the next Section. 
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and 



P!:{v, direct) ^ Ph PI'(v) = P^ Pl^{uh^2 ^ P'/A^. mverse)i^2 , (68) 

as previously observed in Fig. 5. We remark that the above Eqs. ( |67D and ( |68| ) represent 
approximate mirror symmetries, which become exact only for f/^g = 0. However, for U^^ 
obeying the bound in Eq. (|37D , such symmetries are broken, at most, at the few percent 
level. 



IV. INCLUDING EARTH MATTER EFFECTS 

We briefly review known analytical results about Earth matter effects, for the sake of 
completeness and self-consistency of the paper. For recent phenomenological studies of such 
effects in the supernova neutrino context see, e.g., 0-^|l8|-|2^. 

In general, possible Earth matter effects preceding supernova v detection can be imple- 
mented by the final-state substitution [Q 

{Ul,, U!„ Ul) ^ (Pel, Pe2, Pes) (69) 

in Eq. (|60|) or (|65|), where P^i = P{ve ^i) in the Earth. 



Under the assumption of mass spectrum hierarchy (either direct or inverse), the 3z/ 
calculation of Pei is further reduced to a 2z/ problem 0, 

(Pel, Pe2, Pes) ^ [cOS^ 0(1- P^), COS^ P^, siu^ 0] , (70) 

where 

Pe = Pll (71) 

The task is thus reduced to the 2z/ calculations of PeW) and Pe^P)) which are independent 
on ±m^ and on the hierarchy type (direct or inverse). Analytical expressions for Pe can be 
given for particularly simple (or approximated) situations of Earth matter crossing. 



A. One shell 

If the V trajectory crosses only the Earth mantle, characterized by an approximately 
constant (average) density, Pe{^) is simply given by 

2 \ 2 ( sin2co' \ 

Pe[^) = sin uj + sin2ci;m sin(2c<Jm — 2u;) sin - — — L (72) 

\ 2 sin ZLUm I 



where L is the total pathlength in the mantle, and ujm is defined as in Eqs. (^) and (^31) with 
the appropriate potential V in the mantle [Eqs. (0) and (|^)]. The antineutrino probability 
is obtained through the substitution 

PJV) = Pe(u) , (73) 



16 



leading to 



Psiy) = sin^ uj + sin 2uJm sin(2co'm — 2uj) sin 



2 / sin 2^7 
\ 2 sin 2uJrr, 



L 



(74) 



where oJm is defined through Eqs. (|5ll) and (^). 

Figure 5 shows an example of Earth matter effects on P^^ for mantle crossing (dotted 
curves), assuming p = 4.5 g/cm^, Ye = 0.5, and L = 8500 km (parameters which are 
of interest for SN1987A phenomenology; see, e.g., 0,^]). The neutrino potential in the 
mantle is then Vm = 1-7 x 10~^ eV^/MeV, implying strong effects for Vm ~ and thus 
for (5m^ ~ 2EV ~ O(10~^) eV^, which may be relevant in connection with the so-called 
LMA solution to the solar neutrino problem as widely discussed recently [[f|-|9|JT8|-|20 . 
Notice that the approximate symmetries in Eqs. (|67D and (|68D are preserved by Earth matter 



effects, since they act mainly on the two-neutrino L transition in hierarchical approximation 
[Eq. 



B. Multiple shells 

Neutrino oscillations across two Earth shells with different densities [mantle (M) + core 
(C)] were considered in in the context of supernovae, and intensively studied in P2|J43[] 
on general grounds, with emphasis on interesting interference properties peculiar to layered 
matter (see also |0). Adapting, e.g., the notation of to ours, the expression of Pe{^) 
for a mantle+core-|- mantle path in the Earth {L = + Lc + Lm) reads 

Pe(z/) = sin^ to + Wi {Wi cos 2uj + W3 sin 2iu) , (75) 

where W^i 3 are the first and third component of the vector 

W = 2SMYnM + Scnc , (76) 

having defined Y as 

Y = Cm Cc - [riM ■ n^) Sm Sc , (77) 

the vectors um and nc as 

^M,c = (sin2u;A/,c', 0, -cos2ujm,c) , (78) 



and Cm and Cc as 



( kL sin 2^ \ 
Cm,c = cos . Lm,c (79) 



2 sin2ct;M. 



c 



(and similarly for Sm,c, with cos sin), where um and uc are the effective neutrino mixing 
angles in the mantle and in the core. PeO^) is then obtained from Pe{v) through the 
replacement indicated in Eq. (^). 

In the general case of different shells, not necessarily with constant density in each 
shell, the calculation of Pe can also be performed analytically, through the perturbative 
approach developed in [EH] in the context of solar neutrinos. 
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V. SUMMARY AND CONCLUSIONS 



In the context of two-flavor (anti)neutrino transitions in supernovae, we have described 
a simple and accurate analytical prescription for the calculation of the survival probability 
Pee, based on a double exponential form for the crossing probability, and inspired by the 
condition of maximum violation of adiabaticity. The prescription holds in the whole oscil- 
lation parameter space and for generic supernova density profiles. The analytical approach 
has then been generalized to cover three-flavor transitions with mass spectrum hierarchy 
(either direct or inverse), and to include Earth matter effects. 

The final prescription for P^^ can be summarized as follows: 

1. Assume a supernova potential profile V{x) > 0; 



2. Fix the mixing angles (f) = 613 [obeying Eq. (|37|) ] and u = 612, and calculate the matrix 
elements [/^ through Eqs. ; 



X=XL,H 



3. Fix the "solar" and "atmospheric" squared mass differences, Sm"^ > and > 0, 
respectively [within the phenomenological restrictions in Eqs. (^3D-(^)]; 

4. At a given (anti) neutrino energy E, find the points xl and xh where the potential V 
equals the wavenumbers ki = 6m'^/2E and kn = m'^/2E, respectively: 

V{xl,h) = kL,H ; 

5. Calculate the corresponding density scale factors and as 

r 1 dVjxy 
V{x) dx 

(both > for monotonically decreasing V^); 

6. Assign the ± signs in P^ and by choosing the neutrino type (z/ or 77) and hierarchy 
(direct or inverse) in Table II, and calculate then Piu through 

exp(±27rri;,A;i cos^ cj) — 1 
^ exp(±27rriA;i) — 1 

and 

exp{±2'7irHkH cos^ </>) — 1 
^ exp(±27rr///cj^) — 1 ' 

7. Calculate P^^ as 

p!: = U'e, Pt + vl^ (1 - Pt) P| + ul (1 - P|) ; 



8. Finally, include possible Earth matter effects as reviewed in Sec. IV. 

We think that such analytical prescription may be useful to simplify the calculation (and to 
help the understanding) of supernova neutrino oscillation effects. 
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TABLES 



TABLE L Coefficients Xi to be used in the parametrization of P^^ given in Eq. (|60|), according 
to the four possible combinations of neutrino types {u or v) and mass spectrum hierarchy (direct 
or inverse). The coefficients agree with those derived in Refs. [0,^]. 



Type 


Hierarchy 


Xi 


X2 


X3 


V 


direct 




[1-Pl{v)]Ph 


1-Ph 


V 


direct 


1-Pl{v) 


Pl{v) 





V 


inverse 









V 


inverse 


[l-PLmPu 


PLiy)PH 


1-Ph 


TABLE IL 


Signs assigned to Pf^ 


and P^, to be used 


in the parametrization of P^^ 


given in 



Eq. (|65|), according to the four possible combinations of neutrino types {ly or v) and mass spectrum 
hierarchy (direct or inverse). 



Type 


Hierarchy 




P^ 


V 


direct 


+ 


+ 


v 


direct 






V 


inverse 


+ 




V 


inverse 




+ 
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neutrino potential profile for M = 14.6 Mq 
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Fig. 1. Neutrino potential profiles V{x) considered in this work. Dashed curve: "realistic" 



potential, as graphically reduced from the supernova simulation performed in |29| assuming I4.6M0 
for the ejecta. Solid line: "power-law" potential {V oc x^"^) which best fits the realistic one. 
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neutrino crossing probability Pc (power-law potential) 
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Fig. 2. Two-flavor transitions: neutrino crossing probability Pc{i^) in the parameter space 
{Am?/E, tan^ ^) for the power-law potential profile in Fig. 1. Dotted curves: exact numerical 
calculations. Solid curves: results of the analytical prescription in Eqs. (jl^)-(|l4|). Isolines of Pc 
for antineutrinos (not shown) can be obtained by reflection around the axis tan^e = 1. See the 
text for details. 
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neutrino crossing probability Pc (realistic potential) 
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Fig. 3. As in Fig. 2, but for the realistic potential profile in Fig. 1. 
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2y survival probability, E = 15 MeV 




tan "d- 



Fig. 4. Two-flavor transitions: analytical results for the electron neutrino survival probability 
Pee in the mass-mixing plane (Am^, tan^^), at a representative neutrino energy {E = 15 MeV). 
Solid curves: realistic potential. Dotted curves: power-law potential. Isolines of Pee for antineu- 
trinos (not shown) can be obtained by reflection around the axis tan^ = 1. See the text for 
details. 
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5iy survival probability (power-law + Earth mantle) 
E = 1 5 MeV, m' = 3 X 1 0'' eV', tonV = 4 x 1 0"' 




10^10'*10"'l0^10^ 1 10 10^ 10^10^ 1 10 10^ lO"' lo"* 10^ 



tan^w tan^w 



Fig. 5. Three- flavor transitions: analytical results for Pee in the mass-mixing subspace 
{5m?, tan^ cj), assuming the power-law profile in Fig. 1 and fixing E = lb MeV, m? = 3 x 10~^ eV^, 
and tan^ (f) = 4:X 10~^ {P^ = 0.48). Dotted lines include Earth matter effects for a representative 
path of 8500 km in the mantle (assuming p = 4.5 g/cm^ and 1^ = 1/2). The upper panels refer to 
V (left) and V (right) in the case of direct hierarchy. The lower panels refer to v (left) and V (right) 
in the case of inverse hierarchy. 
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